Dynamic phase transition properties and hysteretic behavior of a ferrimagnetic 
core-shell nanoparticle in the presence of a time dependent magnetic field 
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We have presented dynamic phase transition features and stationary-state behavior of a ferri- 
magnetic small nanoparticle system with a core-shell structure. By means of detailed Monte Carlo 
simulations, a complete picture of the phase diagrams and magnetization profiles have been pre- 
sented. According to Neel nomenclature, the magnetization curves of the particle have been found 
to obey N-type, P-type and W-type classification schemes under certain conditions. W-type be- 
havior observed in the present work can be regarded as a novel feature of the dynamic behavior of 
the magnetic nanoparticle system, since the equilibrium counterpart of the model does not exhibit 
such kind of phenomenon. Much effort has been devoted to investigation of hysteretic response of 
the particle and we observed the existence of triple hysteresis loop behavior which originates from 
the existence of a weak ferromagnetic core coupling J c /J 3 h, as well as a strong antiferromagnetic 
interface exchange interaction Ji n t/J s h- Particular attention has also been paid on the influence of 
the particle size on the thermal and magnetic properties, as well as magneto-optical features such 
as coercivity, remanence and compensation temperature of the particle. Furthermore, we have not 
found any evidence of the first order phase transitions in the present system. 

PACS numbers: 64.60.Ht,75.10.Hk,75.70.Rf, 75.75.-c 
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I. INTRODUCTION 



In recent years, influences of small-size and surface ef- 
fects on the magnetic properties of magnetic nanoparti- 
cles have provided a conspicuous and productive field for 
the interaction between theoretical works [l[ and tech- 
nological @, as well as biomedical applications 0, 0]. 
As the physical size of a magnetic system reduces to a 
characteristic length, surface effects become dominant on 
the system, hence, some unusual and interesting mag- 
netic phenomena can be observed, which may differ from 
those of bulk materials IS] . Recent developments in the 
experimental techniques allow the scientists to fabricate 
such kinds of fine nanoscaled materials [3, 0], and the 
magnetization of certain nanomaterials such as 7-Fe203 
nanoparticles has been experimentally measured Q. In 
particular, magnetic nanowires and nanotubes such as 
ZnO @, FePt and Fe3C>4 [13] can be synthesized by var- 
ious experimental techniques and they have many appli- 
cations in nanotechnology [ill, EH > and they are also uti- 
lized as raw materials in fabrication of ultra-high density 
magnetic recording media [l3|-[l5j]. 



From the theoretical point of view, many studies have 
been performed regarding the magnetic properties of 
nanoparticles, such as nanowire, nanotube, and nanorod 
systems, and theoretical works can be classified in two 
basic categories based on the investigation of equilib- 
rium or noncquilibrium phase transition properties of 
such nanoscaled magnetic structures. Namely, in the 
former group, equilibrium properties of these systems 
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have been investigated by a varietyof techniques such 
as mean field theoryJMFT) QJ-Ol [UHll, effective- 
field theory (EFT) fia434j |. Green functions formalism 
[3o| . variational cumulant expansion (VCE) [§3, H3] and 
Monte Carlo (MC) simulations (38l - [53j . Based on MC 
simulations, particular attention has been paid on the 
exchange bias (EB) effect in magnetic core-shell nanopar- 
ticles where the hysteresis loop exhibits a shift below the 
Neel temperature of the antiferromagnetic shell due to 
the exchange coupling on the interface region of ferro- 
magnetic core and antiferromagnetic shell. The readers 
may refer to (45| for a rigorous review about the EB phe- 
nomena. 

It is a well known fact that physical properties of a 
bulk material are independent from size; however, below 
a critical size, nanoparticles often exhibit size-dependent 
properties, and some unique phenomena have been re- 
ported, such as superparamagnetism [5o1. |56|. quantum 
tunneling of the magnetization [57| , and unusual large co- 
ercivities [58| . As an example, it has been experimentally 
shown that Lao.67Cao.33Mn03 (LCMN) nanoparticle ex- 
hibits a negative core-shell coupling, although the bulk 
LCMN is a ferromagnet [13, [§3] ■ Moreover, as a theoret- 
ical example, the total magnetizations in a nanoscaled 
transverse Ising thin film with thickness L are inves- 
tigated by the use of both the EFT with correlations 
and MFT, and it has been shown that the magnetization 
may exhibit two compensation points with the increasing 
film thickness [24j . The phenomena of two compensation 
points observed in the nanoscaled thin films has also been 
reported for bulk ferrimagnetic materials [6l| - l65j . How- 
ever, the origin of the existence of such a phenomenon 
in the nanoscaled magnets are quite different from those 
observed in the bulk ferrimagnetic materials. Namely, a 
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compensation point originates in the bulk systems due to 
the different temperature dependence of the atomic mo- 
ments of the sublattices [66]. However, nanoscaled mag- 
netic particles such as nanowires or nanotubes exhibit a 
compensation point, due to the presence of an antifcr- 
romagnetic interface coupling between the core and the 
shell, even if the lattice sites in the particle core and 
shell are occupied by identical atomic moments. Hence, 
theoretical investigation of ferrimagnetism in nanopar- 
ticle systems has opened a new field in the research of 
the critical phenomena in nanoscaled magnetic particles 
[l7j . Compensation point and coercivity of ferrimagnetic 
materials or nanoscaled magnetic particles are of partic- 
ular interest, since these physical quantities are key fea- 
tures and they hav e gr eat importance in magneto-optical 
recording process [67j. According to recent MC studies 
[50l - l53j , it has been shown that the core-shell morphology 
can be successfully applied in equilibrium properties of 
nanoparticles formed by more than one compound (i.e. 
ferrimagnetic nanostructures) since the concept is ca- 
pable of explaining various characteristic behaviors ob- 
served in nanoparticle magnetism. Namely, we learned 
from these works that compensation point fairly depends 
on the particle size. Therefore, nanoscaled magnets such 
as nanowires, nanotubes, etc. are currently considered as 
promising candidates due to their potential utilization as 
ultra-high density recording media. 

On the other hand, a magnetic system exhibits 
nonequilibrium phase transition properties in the pres- 
ence of a driving magnetic field. Namely, when a mag- 
netic material is subject to a periodically varying time 
dependent magnetic field, the system may not respond to 
the external magnetic field instantaneously which causes 
interesting behaviors due to the competing time scales of 
the relaxation behavior of the system and periodic ex- 
ternal magnetic field. At high temperatures and for the 
high amplitudes of the periodic magnetic field, the sys- 
tem is able to follow the external field with some delay 
while this is not the case for low temperatures and small 
magnetic field amplitudes. This spontaneous symmetry 
breaking indicates the presence of a dynamic phase tran- 
sition (DPT) [H-[70| which shows itself in the dynamic 
order parameter (DOP) which is defined as the time av- 
erage of the magnetization over a full period of the oscil- 
lating field. Related to this nonequilibrium phenomena, 
in recent years, based on Glauber type of stochastic dy- 
namics [71|, a few theoretical studies have been devoted 
to the investigation of dynamical aspects of phase tran- 
sition properties of cylindrical Ising nanowire and nan- 
otube systems in the presence of a time-dependent mag- 
netic field within the EFT with correlations (zi,!^. In 
those studies, the authors analyzed the temperature de- 
pendencies of the dynamic magnetization, hysteresis loop 
area and dynamic correlation between time dependent 
magnetization and magnetic field, and it has been re- 
ported that dynamic magnetization curves can be classi- 
fied into well known categories, according to Neel theory 
of ferrimagnetism [74l l75j ]. Furthermore, based on MC 



simulations and by using unixally anisotropic Heisenberg 
model, frequency dispersion of dynamic hysteresis in a 
core-shell magnetic nanoparticle system has been stud- 
ied by Wu et al. (76|, in order to determine whether the 
dynamic hysteresis loops obey the power-law scaling or 
not, and they concluded that the frequency dispersion 
of the dynamic hysteresis shows both the spin-reversal 
and spin-tilting resonances, and also they found that the 
exchange coupling on the core-shell interface has no ef- 
fect on the power-law scaling of the dynamic hysteresis 
dispersion. 

As can be seen in the previously published works men- 
tioned above, equilibrium phase transition properties of 
nanoparticle systems have been almost completely under- 
stood, whereas nonequilibrium counterparts needs par- 
ticular attention and the following questions need to be 
answered: (i) What is the effect of the amplitude and fre- 
quency of the oscillating magnetic field on the dynamic 
phase transition properties (i.e. critical and compen- 
sation temperatures) of the nanoparticle systems? (ii) 
What kind of physical relationships exist between the 
magneto-optical properties (compensation point and co- 
ercivity) of the particle and the system size? Main mo- 
tivation of the present paper is to attempt to clarify the 
physical facts underlying these questions. For this pur- 
pose, the paper is organized as follows: In Section [TTI we 
briefly present our model. The results and discussions are 
presented in Section IIII1 and finally Section IIVI contains 
our conclusions. 



II. FORMULATION 

We consider a cubic ferrimagnetic nanoparticle com- 
posed of a spin-3/2 ferromagnetic core which is sur- 
rounded by a spin-1 ferromagnetic shell layer. At the 
interface, we define an antiferromagnetic interaction be- 
tween core and shell spins (figure [I}. Construction of 
such kind of model allows us to simulate a ferrimagnetic 
small particle formed by more than one compound. The 
particle is subjected to a periodically oscillating mag- 
netic field. The time dependent Hamiltonian describing 
our model of magnetic system can be written as 

H = -Jmt ^ &iSk — Jc ^ CTjCTj - J s h SkSl 

<ik> <ij> <kl> 

where a = ±3/2, ±1/2 and S = ±1,0 are spin variables 
in the core and shell sublattices. Ji n t, J c and J s h define 
antiferromagnetic interface and ferromagnetic core and 
shell exchange interactions, respectively. h(t) — hsin(u>t) 
represents the oscillating magnetic field, where h and u) 
are the amplitude and the angular frequency of the ap- 
plied field, respectively. Period of the oscillating mag- 
netic field is given by r = 2it/uj. (...) denotes the nearest 
neighbor interactions on the lattice. We fixed the value 
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FIG. 1: Schematic representation of a single core-shell 
nanoparticle on a simple cubic lattice with a spin-3/2 core 
and spin-1 shell structure. 



of J s h throughout the simulations, and we also scaled 
the exchange interactions with J s h. Accordingly, ampli- 
tude of the oscillating magnetic field has been scaled as 
ho = h/J s h in the calculations. 

In order to simulate the system, we employ Metropo- 
lis MC simulation algorithm [79| to equation (QJ on an 
L x L x L cubic lattice with free boundary conditions 
(FBC) which is an appropriate choice for such a finite 
small system. Configurations were generated by select- 
ing the sites in sequence through the lattice and making 
single-spin-flip attempts, which were accepted or rejected 
according to the Metropolis algorithm, and L 3 sites are 
visited at each time step (a time step is defined as a MC 
step per site or simply MCS). The frequency / of the 
oscillating magnetic field is defined in terms of MCS in 
such a way that 



/ = 



1 



k6 s ' 



(2) 



where k is the number of MCSs necessary for one com- 
plete cycle of the oscillating field and Q s is the time in- 
terval. In our simulations, we choose S = 1. hence 
we get t = k. Data were generated over 50 — 100 in- 
dependent sample realizations by running most of the 
simulations for 25000 Monte Carlo steps per site after 
discarding the first 5000 steps. This amount of transient 
steps is found to be sufficient for thermalization for al- 
most the whole range of the parameter sets. However, 
for evaluating the hysteresis loops, in order to guarantee 
to obtain stable loops, the first few cycles of the exter- 
nal field are considered as transient regime, and after 
this transient regime, statistical averaging has been per- 
formed (see Section [III B|) . Error bars were calculated by 
using the jackknife method (80j. 

Our program calculates the instantaneous values of the 
core and shell layer magnetizations M c and tf,/, , and the 



total magnetization Mt at the time t. These quantities 
are defined as 



M c {t) 



M T (t) 



1 iV c - N ah 

-5>, M sh (t) = — j>, 



2 — 1 

M c {t)+M sh (t) 



(3) 



where N c and N s h denote the number of spins in core 
and shell layers, respectively. From the instantaneous 
magnetizations, we obtain the dynamic order parameters 
as follows 

Qc = ^<j)M c (t)dt, Q sh = I M sh (t)dt, 

Q = ^jM T {t)dt, (4) 

where Q c , Q s h and Q denote the dynamic order param- 
eters corresponding to the core and shell layers, and to 
the overall lattice, respectively. We also calculate the 
time average of the total energy of the particle including 
both cooperative and field parts over a full cycle of the 
magnetic field as follows [81| 



Efnt — — 



2ttL 3 



<ik> 



<%J> 



+j sh J2 s kSi + Kt) ( J2 <n + J2 Sk ) 

<kl> \ i k ) . 

Thus, the specific heat of the system is defined as 

dE to t 



dT ' 



dt. 

(5) 

(6) 



where T represents the temperature. To determine the 
dynamic compensation temperature T comp from the com- 
puted magnetization data, the intersection point of the 
absolute values of the dynamic core and shell magnetiza- 
tions were found using 



|0c(T ( 



comp ) 



\Qsh(T c 



comp ) | j 



sign(Q c (T comp )) = -sign((5 s/l (T comp )) 



(7) 



(8) 



with T comp < T c , where T c is the critical temperature i.e. 
Neel temperature. Equations ||7J) and ([8]) indicate that 
the signs of the dynamic core and shell magnetizations 
are different, however, absolute values of them are equal 
to each other at the compensation point. 



III. RESULTS AND DISCUSSION 

In this section, we will focus our attention on the dy- 
namic phase transition properties of the ferrimagnetic 
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nanoparticle system. This section is divided into three 
parts as follows: In Section fill A[ we have examined the 
dependence of the critical temperature T c and compen- 
sation point T comp of the particle on the amplitude and 
frequency of the oscillating magnetic field, as well as the 
exchange couplings defined in equation ([T]). Hysteretic 
response of the particle to the periodically oscillating 
magnetic fields have been investigated in Section IIIIB1 
and size dependence of the magneto-optical properties 
have been analyzed in Section IIII CI In order to make 
a comparison with the previously published works where 
the equilibrium properties of the present system were dis- 
cussed, we select the number of core and shell spins as 
N c = ll 3 and N s h — 15 3 — ll 3 (the same values as in 
53]) in Sections [nTA] and HUH 



A. Dynamic phase transition features of the 
particle 

In order to clarify the influence of antiferromagnetic in- 
terface coupling Jint/Jsh between core and shell layers on 
the dynamic phase transition properties of the system, we 
represent the phase diagrams in a (T Cl T comp — Jint/Jsh) 
plane where hollow (solid) symbols denote the compen- 
sation points T comp (critical temperatures T c ) with three 
oscillation period values r = 50, 100, 200, and for some 
selected values of the field amplitude ho in figures [2M2F- 
Here, we consider a weak ferromagnetic interaction, such 
as Jc/Jsh — 0.25 for the core spins. One of the com- 
mon findings in these figures is that both T c and T comp 
values gradually increase as the strength of the antifer- 
romagnetic interface coupling Jint/Jsh increases, except 
Tcomp curves in figures dp, HH, and[2f. This phenomenon 
is independent from amplitude ho and period r of the 
oscillating magnetic field. On the other hand, as an- 
other conspicuous observation, we can clearly see that 
the system exhibits two compensation points in the high 
amplitude region of {T c ,T comp - Jint/Jsh) plane. T comp 
value which is observed at high temperatures exhibits 
a decreasing trend with increasing \Ji n t\/Jsh values for 
r = 50, it becomes almost stable for r = 100, and dis- 
plays a raising behavior for r = 200. These high temper- 
ature T comp curves disappear at a critical J* nt /J s h value. 
At high oscillation period values (i.e. at low frequencies), 
dynamic magnetization Mt corresponding to the instan- 
taneous ferrimagnetic order parameter of the particle can 
respond to the oscillating magnetic field with some de- 
lay whereas as the period of the external magnetic field 
gets lower, a competition occurs between the period r 
of the field and the relaxation time of the system, hence 
the dynamic magnetization cannot respond to the ex- 
ternal field due to the increasing phase lag between the 
field and the magnetization Mt- As a result, this makes 
the occurrence of the dynamic phase transition difficult. 
In addition, for weak \ Ji n t\ values, core and shell layers 
of the particle become independent of each other. As 
the strength of the antiferromagnetic interface interac- 



tion gets increased then it becomes dominant against the 
periodic local fields, and the particle exhibits a strong 
ferrimagnetic order. Hence, a relatively large amount of 
thermal energy is needed to observe a dynamic phase 
transition in the system, due to the response of the spins 
to the external magnetic field. As the value of the field 
amplitude increases then the antiferromagnetic exchange 
interaction Ji n t loses its dominance against the exter- 
nal field amplitude and it becomes possible to observe a 
dynamic phase transition at lower temperatures. Con- 
sequently, ferrimagnetically ordered phase region in the 
phase diagrams shown in figures [^i-[2f gets narrower with 
increasing ho and r values. In figures [2g and[2)i, we de- 
pict the effect of the antiferromagnetic interface coupling 
Jint/Jsh on the temperature dependencies of dynamic 
order parameters, corresponding to the phase diagrams 
shown in figures [2J: and [2]i. We can see from figure [2g 
that the system may exhibit two T comp values depend- 
ing on the value of J;„ t / J s t ■ According to figure [2ji, 
Tcomp value which is observed at low temperatures orig- 
inates from the temperature dependence of the particle 
core magnetization whereas the other one which is ob- 
served at high temperatures is a result of the thermal 
variation of the shell layer magnetization. As an inter- 
esting observation, although the ferromagnetic exchange 
coupling of the particle core is relatively weaker than 
that of the shell layer (i.e. J c /J s h = 0.25), both the 
core and shell layers undergo a dynamic phase transition 
at the same critical temperature which is a result of the 
relatively strong interface coupling Ji n t- As seen in the 
magnetization curves shown in figure [5g, magnetization 
of the present nanoparticle system can exhibit similar 
features as observed in the bulk ferrimagnetic systems. 
In the bulk ferrimagnetism of Neel [ll, El , it is possible 
to classify the thermal variation of the total magneti- 
zation curves in some certain categories. According to 
this nomenclature, the system exhibits N-type behavior 
for Jint/Jsh = —1-5 and —1.75 which is characterized 
by one compensation point at which resultant magneti- 
zation disappears due to the complete cancelation of the 
magnetization of core and shell layers. The curves corre- 
sponding to Jmt/Jsh = —2.0, —2.5, and —3.0 values show 
W-type behavior which exhibit two compensation points 
before the transition temperature. At this point, we note 
that W-type behavior observed in the present work can 
be regarded as a novel feature of the dynamic behavior of 
the magnetic nanoparticle system, since the equilibrium 
counterpart of the model does not exhibit such kind of 
phenomenon [53| . 

Next, in figure [31 we present the influence of the ampli- 
tude ho of the external field on the temperature depen- 
dencies of magnetization, as well as the dynamic heat ca- 
pacity curves of the particle, corresponding to the phase 
diagrams depicted in figures^ and [2b. In figure [3k., dy- 
namic order parameter curves of the overall system are 
plotted. As seen in figure^, both T c and T comp values 
decrease with increasing ho values. It is clear from figure 
[3Jd that compensation point observed in figure [5£i origi- 
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FIG. 2: Phase diagrams of the system in (T c ,T comp — Jint/Jsh) plane for J c /J s h = 0.25 with some selected values of external 
field amplitude ho = 0.0,0.25,0.50,0.75,1.0. The curves are plotted for three values of oscillation period: (a)-(b) r = 50, 
(c)-(d) r = 100 and (e)-(f) r = 200. Hollow and solid symbols denote compensation point T comp , and critical temperature T c , 
respectively, (g) and (h) show the effect of the antiferromagnetic interface coupling between core and shell spins on the total, 
core and shell magnetizations with temperature for a combination of Hamiltonian parameters. 



nates due to the thermal variation of the magnetization 
of the particle core. On the other hand, dynamic heat 
capacity curves which are depicted in figure [3J; exhibit a 
hump at lower temperatures and a sharp peak at the 
transition temperature. Schottky-like rounded humps 
observed in the heat capacity curves get sharper with 
increasing ho which is a result of a sudden change in 
the core magnetization. Behavior of the dynamic spe- 
cific heat curves corresponding to N-type, and W-type 
magnetization profiles can be better observed in figure [4] 
Namely, when the magnetization curves exhibit N-type 
behavior, dynamic specific heat curves exhibit a hump, 
and a sharp peak whereas if the magnetization shows W- 
type behavior then the specific heat curves exhibit two 
distinct Schottky-like rounded humps, and sharp transi- 
tion peak. When the amplitude ho of the external field 
is sufficiently large, then the first hump observed at low 
temperatures gets sharper whereas the shape of the other 
one which is observed at higher temperatures does not 
change, since it originates from the thermal variation of 
the magnetization of the shell layer which does not ex- 



hibit a sudden variation as the temperature is varied. 

In figure O we examine the effect of the external field 
period r on the dynamical phase transition characteris- 
tics of the particle. Phase diagrams in this figure are 
plotted for a large value of the field amplitude, since we 
can also observe the behavior of additional T comp curve. 
From figures andO), one can clearly observe that both 
Tcomp and T c values are depressed with increasing r. The 
physical facts underlying the behaviors observed in figure 
[5] are identical to those emphasized in figure [3] Therefore 
we will not discuss these interpretations here. However, 
as a complementary investigation, let us represent cer- 
tain magnetization profiles corresponding to the phase 
diagrams given in figure For instance, as seen in fig- 
ures^ and[5J;, when the antiferromagnetic exchange in- 
teraction strength is selected as Jint/Jsh = —0.5 and 
—2.5, observed magnetization characteristics are N-type 
and W-type, respectively. 

As a final investigation of this subsection, we will dis- 
cuss the influence of the ferromagnetic exchange interac- 
tion of the core layer J c / J s h on the magnetization profiles 
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T/J sh 

FIG. 3: Effect of the external field amplitude ho on the tem- 
perature dependencies of (a) total magnetization Q, (b) core 
and shell magnetizations Q c , Q s h and (c) dynamic heat ca- 
pacity C for Jint/Jsh = —0.5, Jc/Jsh = 0.25 and r = 50 with 
h = 0.0,0.25,0.50,0.75,1.0. 
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FIG. 5: Phase diagrams of the particle in (T c ,T cornp — 
Jint/Jsh) plane for J c /J s h = 0.25, ho = 1.0 and r = 
50, 100, 200. Hollow and solid symbols denote compensation 
point Tcomp, and critical temperature T c , respectively. 
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FIG. 4: Temperature dependencies of total magnetization Q, 
core and shell magnetizations Q c , Qsh, and specific heat C of 
the particle for J c /J s h = 0.25, ho = 1.0 and r = 50 with (a) 
Jint/Jsh = -0.5 and (b) Jint/Jsh = -2.5. 



FIG. 6: Effect of the oscillation period r on the temperature 
dependencies of total, core and shell magnetizations of the 
particle for J c /J sh = 0.25, h = 1.0 and r = 50, 100, 200 with 
Jint/Jsh = -0.5 in (a), (b), and Jint/Jsh = -2.5 in (c), (d). 
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FIG. 7: Effect of J c /Jsh on the temperature dependencies 
of total, core and shell magnetizations of the particle for 
Jint/Jsh = -1.0, r = 50 and ft = 0.25,0.5,0.75,1.0 with 
Jc/Jsh = 0.3 in (a), (b), and J c /J s h = 0.6 in (c), (d). 
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FIG. 8: Variation of the coercive field as a function of the 
oscillation period r for several values of ho with some selected 
values of exchange interaction couplings Jint/Jsh and J c /J s h- 



of the particle in figure [7J At first sight, according to fig- 
ure[71 we can clearly claim that as the value of J c /J s h in- 
creases then the transition temperature of the system also 
increases. We can also mention that the phase transition 
temperature of the particle shell is directly related to the 
value of Jc/Jsh, since the antifcrromagnetic interface in- 
teraction is relatively large as Jint/Jsh = —1-0. Further- 
more, for Jc/Jsh = 0.3, the magnetization of the particle 
exhibits N-type behavior, whereas for J c / J s h = 0.6, com- 
pensation point disappears, and we find a P-type depen- 
dence at which the magnetization shows a temperature- 
induced maximum with increasing temperature. These 
observations are represented in figures [7^, and |7£, respec- 
tively. 



B. Hysteretic response of the particle to the 
periodically oscillating magnetic fields 

Hysteresis behavior in magnetic systems originates in 
response to varying magnetic fields, and it is one of the 
most important features of real magnetic materials. In 
dynamic systems, the phenomenon occurs as a result of a 



dynamic phase lag between instantaneous magnetization 
and periodic external magnetic field. In contrast to the 
behavior observed in static models where the strength of 
the external field does not change with time explicitly, 
dynamic hysteresis in nonequilibrium phase transitions 
is characterized by a dynamic symmetry loss at high os- 
cillation frequencies of the external field. The shape of a 
hysteresis loop is determined by the coercivity and rema- 
nent magnetization of the magnetic material. In particu- 
lar, coercivity -which is defined as the required amount of 
the external magnetic field to reduce the magnetization 
of a material to zero- is an essential physical property of 
magnetic materials which has a significant importance in 
magneto-optical recording. Moreover, it is worth to note 
that hysteresis loops of equilibrium systems exhibit coer- 
civity in ferromagnetic phase [HiJ whereas coercive fields 
in nonequilibrium systems driven by an oscillating field 
are always observed in the dynamic paramagnetic phase. 

Oscillation period and external field amplitude depen- 
dence of the coercivity of the particle and dynamic hys- 
teresis loops are shown in figures [51II2I which have been 
calculated at a temperature T = 0.8T® where Tj? is the 
transition temperature of the system in the absence of 
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the external field. This choice of the temperature allows 
the system to undergo a purely mechanical phase tran- 
sition (i.e. magnetic field induced transition). In order 
to acquire a stationary behavior, the first 100 cycles of 
the external field has been discarded and the data were 
collected for 400 cycles. 

From figures [8^ and [HJd, for weak ferromagnetic core 
coupling values such as J c / J s h — 0.25, coercivity curves 
exhibit sudden variation with increasing r whereas for 
sufficiently high r they exhibit a stable profile. By com- 
paring figures [8^l and [SJd we see that as the strength 
of antiferromagnetic interface interaction Jint/Jsh in- 
creases then the width of the hysteresis loops get nar- 
rower. Moreover, it clear from figures |5k ar >d [B}d that 
at low oscillation periods, higher amplitude values mean 
large coercive fields, however at high t values, coerciv- 
ity becomes independent from h . However, according 
to figure \Ejp, as the strength of ferromagnetic exchange 
interaction J c /J s h between the core spins increases, co- 
ercivity may reduce to zero at high r values for relatively 
small ho values. On the other hand, as seen in figure [5Ji, 
under certain conditions, oscillation period dependence 
of the coercivity may exhibit some unusual behavior. 
Namely, for the curves corresponding to Jmt/Jsh = —2.5, 
Jc/Jsh — 0.5, and h = 1.0 and 2.5, as r increases then 
the coercive field of the system tends to decrease up to a 
certain r value. As the oscillation period of the external 
field increases further, coercivity of the particle increases. 

In the following analysis, let us investigate some typi- 
cal examples of the dynamic hysteresis loops correspond- 
ing to figure [H] In figure |H1 we represent the hys- 
teresis loops for some selected values of oscillation pe- 
riod, and for ho — 1.0 corresponding to figure [5^ where 
there exists a relatively weak interface interaction such 
as Jmt/Jsh = —0.5 between core and shell layers. From 
figure [9l we see that remanent magnetizations of the 
core and shell layers do not change as t varies. On 
the other hand, width of the loops becomes narrower 
but does not vanish for high r values. By comparing 
figures IH1 and [TU] we observe that although increasing 
the strength of Jmt/Jsh interactions between core and 
shell spins enhances the coercivity of the core and shell 
magnetizations, the total magnetization curves exhibit 
lower coercive fields. Corresponding to ho = 1.0 curve 
in figure [S}d, for relatively strong antiferromagnetic in- 
terface couplings such as Jmt/Jsh = —2.5 the system 
may exhibit another interesting phenomenon which can 
not be predicted from figure [5}d. Namely, as shown in 
figure 110b . the system exhibits triple hysteresis loops 
with an apparently wide middle loop for relatively strong 
Jint/Jc values. We note that observation of such hys- 
teresis behavior is possible at sufficiently low oscillation 
frequency values. More clearly, we have not observed 
this type of behavior in the system at oscillation periods 
t < 200 with the parameters given in figure [TUb . Ori- 
gin of this phenomena can be understood by analyzing 
the time series of instantaneous magnetizations. Accord- 
ing to our numerical calculations, triple hysteresis loop 



behavior originates from the existence of a weak ferro- 
magnetic core coupling J c / J s h, as well as a strong an- 
tiferromagnetic interface exchange interaction Jmt/Jsh- 
Triple loops disappear for Jmt/Jsh > 0. In figure HT1 
hysteretic response of the instantaneous magnetization 
curves to the periodic external field is presented where 
the situation is depicted for when the ferromagnetic core 
coupling of the system is relatively weaker than that of 
the shell layer, i.e. J c /J s h — 0.25, and in the presence 
of both weak and strong antiferromagnetic interface in- 
teractions, i.e. Jmt/Jsh = -0.5 and J lnt /Jsh = -2.5, 
respectively. Oscillation period is selected as r = 10000, 
hence the instantaneous magnetizations are able to fol- 
low the external field almost simultaneously. According 
to figure [TTk where we consider a weak interface coupling 
(Jint/Jsh = —0.5); at strong fields, such as for ho = — 1.0 
both the core and shell layers are magnetized along the 
field direction. As the field amplitude reaches to the value 
ho = 1.0 within a half cycle, both core and shell mag- 
netizations easily keep their alignment with the external 
magnetic field. On the other hand, when Ji„ t / J s h — — 2.5 
(see figurelTTb) , shell magnetization can align in the oscil- 
lating field direction instantaneously, however core mag- 
netization tends to align antiferromagnetically, due to the 
existence of a strong Ji nt interaction and weak ferromag- 
netic core interaction (J c / Jsh — 0.25). This results in a 
maximum phase lag between the external magnetic field 
and magnetization of the core layer, and consequently we 
observe triple hysteresis loops. These types of hysteresis 
loops have also been observed recently in cylindrical Ising 
nanowire systems [28| . and in molecular based magnetic 
materials [83| in the presence of static magnetic fields. 

In figure 1121 we represent hysteresis loop profiles for 
various oscillation periods corresponding to phase dia- 
grams in (coercivity — r) plane in figures |8j; and [8jl. 
Here we consider a high magnetic field amplitude, such 
as ho = 5.0 with J c /J s h — 0.5. As shown in figure [T2"b. 
when we select Jmt/Jsh — —0.5, hysteresis loops get nar- 
rower with increasing r values, however the loop areas 
do not disappear, since the coercive fields do not vanish, 
even at large r values (see figure [5fc). On the other hand, 
from figure W2b , we observe that when the strength of 
antiferromagnetic interface interactions increases, i.e. for 
Jint/Jsh = —2.5, loop areas tend to disappear with in- 
creasing t values, since the coercive field value decreases, 
and reduces to zero at a certain r value. 



C. Size dependent properties 

As a final investigation, let us represent some size de- 
pendent properties of the particle for some selected values 
of Hamiltonian and magnetic field parameters. In figure 
Q21 we show the effect of the shell thickness L s h on the 
coercivity and hysteresis curves corresponding to high 
frequency regime (r = 50). It is clear from figure [T5k 
that as the temperature increases then coercivity of the 
system increases and saturates at a certain value which 
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FIG. 9: Hysteresis loops of the particle with a relatively weak Jint/Jsh value corresponding to (a) core, (b) shell and (c) total 
magnetizations for several values of oscillation period r with some selected values of exchange interaction couplings Jint/Jsh 
and Jc/Jsh- The letters accompanying each curve denote the value of r. 




FIG. 10: Hysteresis loops of the particle with a relatively strong Jint/Jsh value corresponding to (a) core, (b) shell and (c) total 
magnetizations for several values of oscillation period r with some selected values of exchange interaction couplings Jint/Jsh 
and Jc/Jsh- The letters accompanying each curve denote the value of r. 
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FIG. 11: Low frequency (r = 10000) hysteresis loops of the 
particle for a selected field amplitude ho = 1.0 with some se- 
lected values of exchange interaction couplings Jint/Jsh and 
Jc/Jsh- (a) Weak, (b) strong antiferromagnetic interface cou- 
pling values. Solid, dashed and dashed-dotted lines corre- 
spond to core, shell and total magnetization curves, respec- 
tively. 
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FIG. 12: Hysteresis loops of the total magnetization of the 
particle for a field amplitude value ho = 5.0 with several 
values of oscillation period r corresponding to some selected 
values of exchange interaction couplings Jint/Jsh and Jc/Jsh- 
(a) Weak, (b) strong antiferromagnetic interface coupling val- 
ues. The letters on the curves represent the value of t. 



depends on the temperature. Moreover, the curve cor- 
responding to T — 2.5 (black full squares) exhibits an 
unusual behavior. Namely, coercivity of the system loses 
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FIG. 13: (a) Dependence of the coercivity of the particle as 
a function of shell thickness L a h with some selected values of 
temperature. The system parameters have been kept fixed 
as Ji n t/Jsh = -0.1, Jc/Jsh = 0.25, r = 50 and h = 0.5. 
Corresponding hysteresis curves of the total magnetization 
Mt have been depicted for temperature values T = 2.5 and 
3.0 in (b) and (c), respectively. 



its symmetric shape with increasing L s h when T = 2.5 
which can be regarded as a signal of a dynamic phase 
transition. Corresponding hysteresis curves have been 
depicted in figures IT3"b andlTSb for T = 2.5 and T = 3.0, 
respectively. As seen from figure H3b . hysteresis loop 
loses its symmetry and the particle exhibits a dynamic 
phase transition from paramagnetic to a dynamically or- 
dered phase with increasing L s h values. This phase tran- 
sition completely originates from the magnetization of 
the shell layer of the particle, since the magnetization of 
particle core does not change its shape as L s h varies. On 
the other hand, according to figure [T5b . the particle al- 
ways remains in the paramagnetic phase as L s h increases, 
since the temperature is large enough to keep the system 
in a dynamically disordered state. We can also observe 
from figure [T3"b that loop areas get wider with increasing 
L s h, an d after a sufficiently large value of L s h, the loop 
areas are not affected from changing L s h values. We have 
also investigated the situation for low frequency regime 
(t = 5000) in figure [JJJ At low frequency values, the 



magnetization of the particle is able to follow the ex- 
ternal magnetic field instantaneously, but with a phase 
lag. Dependence of coercivity on the shell thickness L s h 
for r = 5000 can be seen in figure [Ha . By compar- 
ing figures H3k and [14k . we can clearly observe that the 
system exhibits large coercivities for high frequency val- 
ues. Hysteresis curves corresponding to figure [T4a are 
plotted in figures [Tib and [Tib for temperatures T = 2.5 
and T = 3.5, respectively. The most remarkable observa- 
tion in these figures is that the phase difference between 
the magnetization of the particle and oscillating external 
magnetic field drastically reduces, hence coercivity and 
remanent magnetization values become zero in low fre- 
quency and high temperature regions, and consequently 
loop areas reduce to zero. 
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FIG. 14: (a) Dependence of the coercivity of the particle as 
a function of shell thickness L s h with some selected values of 
temperature. The system parameters have been kept fixed 
as Jint/Jsh = —0.1, Jc/Jsh = 0.25, r = 5000 and ho = 0.5. 
Corresponding hysteresis curves of the total magnetization 
Mt have been depicted for temperature values T = 2.5 and 
3.5 in (b) and (c), respectively. 

In figures [T51T71 in order to depict the effect of the 
ferromagnetic shell layer thickness L s h on T comp and T c 
values, we plot the temperature dependencies of mag- 
netization curves of the particle for some selected val- 
ues of L s h- Figure [T5k shows the representative W-type 
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FIG. 15: Dependence on the shell thickness L s h of thermal 
variation of (a) total magnetization, (b) core and shell mag- 
netizations for Jint/Jsh = —3.0, Jc/Jsh = 0.25, ho — 1.0 and 
r = 100. 



FIG. 16: Shell thickness L s h dependence of (a) total magneti- 
zation, (b) core and shell magnetizations for Ji„t/ J s h = —0.5, 
Jc/Jsh = 0.25, h = 1.0, and r = 50. 
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curves for some selected values of Hamiltonian parame- 
ters. It is clear from figure [15] that lower T comp value is 
independent of L s h, since this compensation point origi- 
nates from thermal variation of particle core magnetiza- 
tion. However, higher T comp value exhibits a size depen- 
dent behavior. The magnetization curves corresponding 
to L s h = 2 prominently undergo a shape variation as L s h 
increases, and the temperature induced maxima of total 
magnetization curves get higher with increasing thick- 
ness. Similarly, in figure [TBI L s h dependence of represen- 
tative N-type magnetization curves are shown. Accord- 
ing to figure ITBa . although T comp value is not affected 
from changing L s h values, transition temperature of the 
particle increases with increasing L s h when the magneti- 
zation exhibits N-type behavior. Finally, the situation is 
depicted for representative P-type curves in figure [171 It 
can be easily observed in figure ITTa that the temperature 
induced maxima of the total magnetization curves get 
higher with increasing L s h- One of the common prop- 
erties observed in figures [T5l[T7l is that the magnetiza- 
tion curves retain their shapes even if the particle size 
changes. Furthermore, in contrast to the results given in 
figure the curves depicted in figures [15] and [T7] ap- 
proach to bulk limit with increasing L s h values. 
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FIG. 17: Effect of the shell thickness L s h on the temperature 
variation of (a) total magnetization, (b) core and shell mag- 
netizations of the particle for Jint/Jsh = —1.0, J c / Jsh = 0.6, 
h = 0.75, and r = 50. 



change interaction in the interfacial region we have also 
investigated some of the ferrimagnetic properties of the 
particle. The most conspicuous observations reported in 
the present paper can be briefly summarized as follows: 



IV. CONCLUDING REMARKS 

In conclusion, by means of MC simulations, in order 
to clarify how the magnetism in a nanoparticle system is 
affected in the presence of a periodically oscillating exter- 
nal perturbation, we have analyzed nonequilibrium phase 
transition properties and stationary-state behavior of a 
single domain ferrimagnetic nanoparticle which is com- 
posed of a ferromagnetic core surrounded by a ferromag- 
netic shell layer. By considering an antiferromagnetic ex- 



• In Section HII A[ a complete picture of the phase di- 
agrams and magnetization profiles have been pre- 
sented. We have found that the system may exhibit 
two compensation points in the high amplitude re- 
gion of (T c ,T comp - Jint/Jsh) plane, and the high 
temperature T comp curves in the same plane disap- 
pear at a critical J* nt /Jsh value. According to Neel 
nomenclature, the magnetization curves of the par- 
ticle are found to obey N-type, P-type and W-type 
classification schemes under certain conditions. W- 
type behavior observed in the present work can be 
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regarded as a novel feature of the dynamic behav- 
ior of the magnetic nanoparticle system, since the 
equilibrium counterpart of the model does not ex- 
hibit such kind of phenomenon (53[. 

• Section MI Bl has been devoted to investigation 
of hysteretic response of the particle. Based on 
the simulation results, we have reported the exis- 
tence of triple hysteresis loop behavior which orig- 
inates from the existence of a weak ferromagnetic 
core coupling J c /J s h> as wen as a strong antiferro- 
magnetic interface exchange interaction JmtjJsh- 
Triple hysteresis loops disappear under the trans- 
formation J int -> -Jtnt- 

• Size dependent properties of the particle have been 
clarified in Section llll Cl and we found that the par- 
ticle may exhibit a dynamic phase transition from 
paramagnetic to a dynamically ordered phase with 
increasing ferromagnetic shell thickness L s h- 

All of the observations outlined above show that the 
shape (amplitude and frequency) of the driving field and 
particle size have an important influence on the thermal 
and magnetic properties, as well as magneto-optical fea- 
tures such as coercivity, remanence and compensation 
temperature of the particle. In fact, we note that ac- 
cording to our simulation results, we have not found any 



evidence of the first order phase transitions. The reason 
is most likely due to the fact that, in contrast to the con- 
ventional techniques such as MFT and EFT, the method 
we used in the present work completely takes into account 
the thermal fluctuations in the present system which al- 
lows us to obtain non-artificial results. Moreover, it is 
possible to improve the proposed model to simulate more 
realistic systems by considering a simulation of Heisen- 
berg type of Hamiltonian with an assembly of interacting 
nanoparticles instead of a single particle. This may be 
the subject of a future work. We hope that the results ob- 
tained in the present paper can be helpful to experimen- 
tal investigations of nanoscaled magnetic systems with a 
core-shell structure. 
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